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Abstract 


Harmonic oscillator wave functions have played an historically important role in our 
understanding of the structure of the nucleon, most notably by providing insight into the 
mass spectra of the low-lying states. High energy scattering experiments are known to give 
us a picture of the nucleon wave function at high-momentum transfer and in a frame in 
which the nucleon is travelling fast. This paper presents a simple model that crosses the 
twin bridges of momentum scale and Lorentz frame that separate the pictures of the nucleon 
wave function provided by the deep inelastic scattering data and by the oscillator model. 


1 Introduction 


While a prediction for the nucleon structure functions from first principles in quantum chromodj- 
namics (QCD) seems, even now, a remote prospect, if we are ever to claim a deep understanding of 
the structure of the nucleon, a clear interpretation of such gross properties as the neutron-proton 
structure function ratio (R np ) and the polarization asymmetry (A™) of the proton structure func- 
tion is essential. A notable attempt to relate these features to the nucleon rest-frame wave function 
was made by Le Yaouanc et al. [1, 2, 3], who employed non- relativistic harmonic oscillator spatial 
wave functions and SU(6) mixing in an attempt to formulate predictions both about the structure 
functions and the nucleon form factors. While both the large-x behavior of R np and the initial 
slope of the neutron electric form factor were well accounted for by the inclusion of an admixt ure 
70 excited state in the nucleon wave function, the signs of the mixing angles obtained in the two 
cases were observed to disagree. 

Against the structure- functions calculation of Le Yaouanc et al. may be raised the objection that 
there is no clear prescription for Lorentz- transforming a non- relativistic wave function. It is tins 
concern that will be addressed in this paper. Less widely recognized is an objection that can 
be raised against the treatment of the form factors by Le Yaouanc et al. The latter calculation 
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involves the assumption that the nucleon’s charge (magnetization) density and electric (magnetic) 
form factor are related by Fourier transformation. The Fourier relationship holds only when the 
Lorentz transformation of the spatial wave functions is ignored. There are several models which are 
known to predict a non-zero neutron electric form factor in the absence of SU(6) mixing[4, 5, 6]. In 
such models, which employ plausible relativistic spin wave functions, the matrix elements involved 
in the determination of the form factors cannot be factorized into a product of spin, isospin, and 
spatial matrix elements. Since the spin wave functions play no role in the structure function 
calculation, the possibility must be considered that the structure functions provide the correct 
value for the mixing angle. 

The spatial wave functions that shall be considered are the “definite metric” [7] solutions to the 
relativistic harmonic oscillator equation of Feynman et al. [8] In their original work, Feynman et al. 
used the non-normalizable “indefinite metric” solutions of their wave equation. These solutions 
yield divergent form factors as —q 2 increases. To remedy this, they multiplied all matrix elements 
by an ad hoc factor. The “definite metric” solutions are normalizable and, when used to calculate 
nucleon form factors, yield the proper q 2 behavior, a dipole fall-off for large — <? 2 , without any 
adjustments. These solutions also help to illuminate features of the structure functions and the 
parton model, as will be seen later on. 

In Section 2 the relativistic harmonic oscillator equation and its normalizable solutions are re- 
viewed. The behavior of these solutions under Lorentz’s transformation is discussed, and their 
form in the infinite momentum frame is exhibited. In Section 3 the infinite- momentum- frame 
relativistic-oscillator nucleon wave function is combined with QCD momentum scaling incorpo- 
rated via the valon model of Hwa.[9] The proton and neutron structure functions are considered 
within the context of the resulting model, and a value for the mixing angle for an admixture of 
70 excited state is calculated. In Section4, the significance of this calculation is reviewed. 


2 The Relativistic Oscillator Model 


For simplicity of discussion, the relativistic oscillator model is introduced for the two particle case 
first. This model describes the binding of a pair of quarks to from a meson via the differential 
equation 


{2 [d 2 + d 2 ] - (w 2 /16) (ii — x 2 ) 2 — mQ j X2) = 0 (1) 

where Xj and x 2 represent the space-time coordinates of the two constituent quarks, and the metric 
convention is defined by -500 = 9i, = 1 • The quark spin will be ignored here, though versions of 
the relativistic oscillator model have been formulated to include spin 1/2 quarks. [10, 11] Eq. (1) 
is readily solved via separation of variables in terms of the coordinates 
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( 2 ) 


^ fi) 

l/2(xi /i X2 11 ) 


where the X ^ are the space-time coordinates of the meson center of momentum and the x M 
determine the space-time separation of the quarks. The separated equations are 


(d\ - m 2 ) ip(X) = 0 (3) 

and 

( ~d 2 x + w 2 /4x 2 + <^(x) = m 2 (p{x) (4) 

where , J , (xi,X 2 ) = ift(X)<p(x). Eq. (3) is the Klein-Gordon equation for a meson of mass m, while 
Eq. (4) describes a four-dimensional harmonic oscillator. 

Eq. (4) is itself separable in terms of the space-time components x M , while the eigenvalue m 2 is 
given by a linear combination of the eigenvalues corresponding to each of the component equations. 
In the timelike direction, an increase in the excitation quantum number corresponds to a more 
negative contribution to the mass squared. To eliminate a degree of freedom which is not observed 
in nature, and to eliminate, as well, the unphysical possibility of imaginary mass, the oscillator- 
model solutions are required to obey the subsidiary condition 


d 

dX, 


_1 d_ 

dx ** 



^(X,x) = 0. 


( 5 ) 


This condition suppresses timelike excitations in the meson rest frame. 

The solutions to Eq. (3) have the familiar form exp (iP^X^) where is the four-momentum 
corresponding to the meson center-of-momentum coordinates X M . The solutions to Eq. (4) are 
products of oscillator solutions in each of the space-time components, with the solution in the 
timelike coordinate in the restframe being restricted to the fundamental mode via Eq. (5). Such 
solutions can be written as 


^•“V) = (2 6+fc+ -6!A:!^!) _1/2 ^4x r vW2]^[^VW2]^[^VT72] 

x exp [-w/4 (x 2 + y 2 r + z 2 + i 2 )] 


( 6 ) 
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where H denotes a Hermite polynomial and x' denotes the four vector x, represented in terms 
of its components, x r , y T , z r and t r in the meson rest frame. The invariant meson square mass 
corresponding to P M is required by Eq. (3) to be equal to m 2 , while Eq. (4) determines that 
m 2 = u>(6 + fc + «' + l) + rnl- 

The above solutions form a complete set of normalized rest-frame solutions. The wave function 
in a frame in which the meson is not at rest is specified by the Lorentz transformation between 
the meson rest frame and the frame in which it is moving. For example, the ground state in an 
arbitrary frame can be written as 


$(X,x) = w/(2jt) (2 h+k+w b\klw^ 1/2 e' 7 ”'*" 

x exp{-w/4[x 2 - 2(P • x) 2 /P 2 ]} . (7) 


The construction of relativistic-oscillator momentum-space wave functions in arbitrary frames is 
equally straightforward. Figure 1 provides a pictorial view of the effect of the Lorentz transfor- 
mation on the rest-frame wave function, both in coordinate space and in momentum space. The 
bound state quarks are seen to acquire lightlike momenta in the frame where the meson is moving 
rapidly. The success of the parton model tells us that this should be the case. 

Modelling of the nucleon requires that a three particle version of the relativistic oscillator be 
considered. A harmonic interaction between each pair of quarks is assumed, and the governing 
differential equation takes the form 


{3 [d 2 + d\ + dl] - o> 2 /36 [(*i - x 2 ) 2 + (x 3 - x 2 f + (x 3 - x,) 2 ] - U 0 } 9(x u x 2 ,x 9 ) = 0 (8) 

where Xi, x 2 and x 3 are the space-time coordinates of three constituent quarks. 

Separation of variables can be implemented in terms of the coordinates X , r and s, defined as 


X ^ — 1/3 (x ltl T x 2m + x 3 „) 

= l/6(x 1#1 T x 2m 2x 3 ^) 

s M = -1/(2 v / 3)(x 1 ^-x 2#1 ). (9) 
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FIG. 1 The Lorentz deformation properties of the relativistic oscillator in coordinate 
space and in momentum space. 
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The wave function ^(xi, x 2 , x 3 ) can be written in terms of these variables as tp(X)<p(r)0(s) where 
ip(X) satisfies Eq. (3) while <f>(r) and 0(s) satisfy respectively 


-l/2(d r 2 - W V) 9 (r) = A r <^ 

-l/2(d 2 ,-u 2 s 2 )d(s) = X 3 6(s). (10) 


The square mass is Eq. (3) in the nucleon case is then given by U° + X r + X 3 . To remove the 
unphysical timelike degree of freedom from the nucleon spectrum of states determined by Eqs. (10), 
the three-particle relativistic oscillator equation is supplemented by a pair of subsidiary conditions 
that suppress such excitations in the nucleon rest frame. 

Application of the relativistic oscillator model to the structure functions requires construction of 
the momentum-space wave function in the infinite momentum frame. In a frame in which a meson 
whose rest-frame wave function is given by (6) is travelling with velocity parameter (3 along the 
^-direction, the internal momentum-space wave function is 


V W (P,P) 



x exp [— l/u;((p* - 0Po) 2 + (p 0 - 0Pzf)/{ 1 - /? 2 )] 


( 11 ) 


where p represents the momentum conjugate to the internal separation coordinate x, and where 
transverse degrees of freedom have been neglected. As /3 — * 1, the square magnitude of pip, 3) 
becomes singular along the forward light cone, while vanishing everywhere else. Integrating along 
the direction perpendicular to the forward light cone results in a distribution for the internal 
light-cone momentum p + (= p 0 + p z ) given by 

p(p+) = lim J dp. | <p w (p,/3) \ 2 (12) 

where p_ = p 0 - p z . The distribution p(p + ) is converted into a distribution in Feynman x by 
setting pi+ = x P and requiring p(x)dx = p(p + )dp + . 

A similar procedure may be followed in the three-particle case. For three particles the result is[12] 


p{x) = 3m/(2TTuj) 1/2 'jr ^ j (l/il)(l/2Y H 2 [(m/ y/&J)(l - 3x)] 

x exp[— (m 2 /2u;)(l - 3x) 2 ] (13) 
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in general, and 


Po{x) = 3m/(27ro;) 1/ ' 2 • exp[-(m 2 /2u;)(l — 3x) 2 ] (14) 

if the nucleon is assumed to be described by the oscillator ground state wave function. The variable 
x in Eqs. (13) and (14) is the momentum fraction variable. A calculation of the proton charge 
structure function F 2 ep (x), can, for example, be based directly on Eq. (14). The result is 


F 2 p (x ) =< e 2 > mx/(27ru;) 1/,2 • exp[(— 9m 2 /2u;)(x — 1/3) 2 ] (15) 

where the average of the charge e, is taken over the three valence quarks. This calculation ignores 
scaling effects predicted by QCD and yields only qualitative agreement with experiment. 


3 Structure Functions 


A valon is a bound state or constituent quark whose internal structure is probed in high energy 
interactions. To be completely general, valons of different spin as well as flavor should be differ- 
entiated. Let G v /n{x) represent a momentum- fraction probability distribution for a valon of type 
v ( v representing spin and flavor) in the nucleon N. A nucleonic structure function ^^(x,^ 2 ) 
is expressed in terms of convolutions of G v /w(x) with corresponding structure functions for the 
valons: 


F n (x, Q 2 ) = J2 £ dx’G v/N (x')F v (x/x',Q 2 ). 


(16) 


The Q 2 dependence of the structure functions appears only in F v (x, Q 2 ). QCD evolution Eq. (13) 
for the moments of the struct ure functions are used to express this dependence. According to 
Eq. (16), the moments of a nucleon structure function are given by a sum of products of moments: 


M N (n,Q 2 ) = '£M v/N (n)M v (n,Q 2 ) 

V 


(17) 


where 


M N ’ v (n,Q 2 ) = I' dxx n - 2 F N ' v (x,Q 2 ) 
Jo 


(18) 
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and 


M v/N {n ) = jf dxx n ~ l G v/N (x). (19) 

The evolution equations are the basis for assuming a form for the moments M v (n,Q 2 ) of the 
structure functions F v {x,Q 2 ). The F v {x,Q 2 ) are understood to be determined by the quark 
distributions within the valon, v, which distributions can be broken up into components^ that 
behave as singlets and as nonsinglets under flavor transformation. The moments JVf (n, Q ) are 
correspondingly expressed in terms of singlet and nonsinglet moments, which are defined to be the 
scaling factors governing the evolution of the moments of such quark distributions in lowest-order, 
twist-2 QCD. The nonsinglet moments are given by 


A/jvs(n,<3 2 ) = exp {-d^ss) 


( 20 ) 


while the singlet moments are 


M,{n,Q 2 ) = 1/2(1 + p n )exp(-d n + s) + 1/2(1 - p n ) exp(— <Ts) (21) 


where 


s = In 


/ym 

VMCS/A 2 ) 


)' 


(22) 


The coefficients and p n come from the renormalization group analysis. [13] The con- 

stant, A, is the usual scaling param eter while Q 0 represents the “starting point ’ of the evolution. 

Since valons of different helicity as well as flavor are to be distinguished, four separate valon 
distributions will be required to characterize the nucleon. The corresponding moments are denoted 
as 


Ui(n) = Muy p (n) = Mo T / n (n) 

£> T (n) = M DT / P (n) = Mt/|/n(ri) 

ffj(n) = Mui/ P {n) = M£q/„(n ) 

>\'n) = M D i/ P { n ) = Mui/ n {n) (23) 
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where the symbol f (U denotes that the valon’s helicity is par allel (antiparallel) to that of the 
nucleon, and where the identification of valon distributions within the neutron with the corre- 
sponding isospin- reversed distributions in the proton follows from charge symmetry. In terms 
of the singlet and nonsinglet moments Eqs. (20) and (21) and the valon moments Eq. (23), the 
moments of the nucleon structure functions F 2p (x,Q 2 ) and F 2n (x,Q 2 ) are given by 


M 2p (n, Q 2 ) = 2/9[2f/(n) + D{n)]M,{n, Q 2 ) + l/9[4tf(n) - D(n)]M NS {n,Q 2 ) (24) 


M 2n (n, Q 2 ) = 2/9[2 U(n) + D{n)]M,{n,Q 2 ) - 2/9 [U{n) - D(n)]M NS {n, Q 2 ) (25) 

with 

U(n) = U,{n) + U0n) 

D(n) = D T (n) + Z)j(n). (26) 

It is easily verified that these equations describe the lowest-order twist-2 QCD evolution of the 
moments of F 2p and F 2n from a starting point at which the nucleon is viewed as consisting of its 
three bound-state quarks. 

Eqs. (24) and (25) were used by Hwa[9] in conjuncti on with experimental moments of F 2p and 
F 2n to obtain fitted values for the parameters Q 0 and A. These equations are first order, and will 
therefore not be accurate for low Q 2 . Ideally we would like to evolve the bound-state momentum- 

space wave function from the energy scale Ql at which the nucleon is describable as a bound state 

of its three constituent quarks (with, perhaps, an oscillator-like momentum distribution), out to 
high Q 2 where the structure functions are observed. The fitted parameter Qg is an approximation 
for Ql in the sense that the lowest-order evolution equations are used. This approximation is a 
key feature of the valon model and is discussed in detail in. [14] The goal of Hwa’s fitting procedure 
was to obtain estimates for the functions G v/N (x). In Figure 2, an “average” valon distribution 
obtained in[9] by neglecting spin and flavor dependence is compared with pa{x) given by Eq. (14) 

Let us now introduce a Id component of SU(6) into the nucleon wave function in the form 


vf = [cos04> 0 | 56 >, +(smd/\/2)((f> a | 70 > 0 +(f>p | 70 >^)] • exp(— iP ■ X). (27) 

The (f>'s represent the spatial wave functions; <f > o is the harmonic oscillator ground state, while <f> a 
and <f>p are taken to b e excited states with total harmonic oscillator quantum number n = 2 and 
zero orbital angular momentum. The subscripts a and 0 refer to the two possible types of mixed 
symmetry which are characterized by the behavior of the (three quark) wave function under 
exchange of the first and second quarks. The form of the excited-state component is uniquely 
determined in the oscillator model. The 7Q. state that involves n = 1 oscillator wave functions 
is disallowed because it is of the wrong parity. No other n = 2 state with the same quantum 


519 



numbers as the nucleon interferes with the ground state to produce the SU(6) breaking effects 
that are observed in the structure functions. The wavefunction Eq. (27) leads to spin-and-flavor 
dependent valon distributions of the form 



FIG. 2 A comparison of Hwa’s “average” valon distribution with po(x) defined by 
the infinite-momentum-frame relativistic-oscillator momentum-space wavefunction. 


Guy P {x) = (3 m/V2iru) [5/6 cos 2 9 + sin 2 9 {5/36A(x) + l/3i'(x)} - 2v^6/9 cos 9 sin 9j(x)] 
x exp [— (m 2 /2ui)(l — 3x) 2 J 

Gui/p( x ) = (3m/v^ru>) [l/6cos 2 0 + sin 2 # {1/36 h(x) + l/3t(x)} + VQ/lScos9 sin0?(x)] 
x exp [-(m 2 /2w)(l - 3x) 2 )j 

Gdi/p{x) = (3m/\/27raJ) [l/3 cos 2 9 + sin 2 9 {l/lSh(x) + 2/3 z(x)} + \/6/9cos^sin^'(x)| 
x exp [— (m 2 /2w)(l — 3x) 2 )J 

Gdi/p{x) — (3m/V2wu;) [2/3cos 2 0 + 1/9 sin 2 9 + 2\/6/9cos0sin0?(x)] 

x exp [— (m 2 /2u;)(l - 3x) 2 )] (28) 

where 

h(x) = 43/16 + m 2 /8w(l — 3x) 2 + m 4 /16u/ 2 (l — 3x) 4 
i(x) = 5/8 + m 2 /8w(l-3x) 2 

j{x) = 1/4 — m 2 /4u;(l — 3x) 2 . (29) 

Moments U(n) and D(n) determined from the above distributions were used in Eqs. (24) and 

(25) to obtain fits for experimental moments[15] of F 2p (x) and F 2n (x) derived from the CHIO 
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muon data[16] and SLAC electron data[17] at Q 2 = 22.5 GeV 2 . A somewhat large value of Q 2 
was chosen to minimize target mass and higher twist effects that may be present in the data. The 
ratios R np (x) and A' ip {x) do not app ear to show any appreciable Q 2 dependence. The extension 
of the tails of the distributions into the unphysical regions x < 0 and x > 1 was ignored for 
purposes of computing the moments. The resulting small deviation from the Adler sum rule does 
not appear to lead to noticeable discrepancies. 




FIG. 3 The moments of the nucleon structure functions vs. n as fitted by Eqs. (24) 
- (26) in conjunction with Eq. (29). Fitted moments at 9 = 0° and at 9 = 23.3° are 
represented by the solid curves and are compared with data from[15]. 


The fitted moments were functions of two parameters - the mixing angle 9 and the scaling variable 
s defined in Eq. (22). minimization was used to determine the best fit. The \ 2 function in 
this case cannot be taken as an absolute indication of the quality of the fit due to the statistical 
interdependence among the moments of F 2p and F 2n . l}i 2 was used, rather, as a relative determi- 
nant of merit, so that the quality of the fit as a function of 9 could be evaluated. The minimum 
of x 2 occurs at 6 = 23.3°, and a positive mixing angle is clearly preferred. Figure 3 compares the 
best-obtainable predicted moments from Eqs. (24) and (25) for 0 = 0° and for 9 = 23.3° with 
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the experimental moments. At 9 = 0°, the fitted moments of F 2n fall outside the error limits for 
large n. With the inclusion of the IQ. state in the wave function at a mixing angle of 23.3 c,rc , 
a simultaneous fit to the moments of F 2p and F 2n appears more reasonable, although the fitted 
moments of F 2n remain somewhat large for large n. 


4 Conclusion 


The simple model presented in this paper falls short of providing us with the ability to draw 
precise numerical correspondences between nucleonic bound state properties and the structure 
function data. The model does, however, address the crucial questions of Lorentz transformation 
and momentum scaling that must be considered if such corespondences are ever to be drawn. The 
approximate agreement between p 0 (z) and Hwa’s phenomenologically-determined valon distribu- 
tion (see Figure 2) allows us to believe that some of the essential physics is being captured. The 
value of the SU(6) 70 state mixing coefhcien t obtained in this model via a simultaneous fit to 
proton and neutron structure function moments is very close to the original value determined by 
Le Yaouanc et ai This fact, together with the dependence of the form factors on the nucleon 
spin wavefunction, lends credence to the idea that the observed behavior of R np and A 1V can be 
reliably interpreted as evidence of SU(6) mixing in the nucleon wavefunction. 
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